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Lunar thermal history model Hng demonstrates that a realistic computer simu~ 
lation can be used to predict many more observables than heat flow and in- 
ternal composition. These observables include major and minor element con- 
centrations, thickness of the lunar "crust", intensity of volcanic activity 
as a function of time, etc. 

The models which are most consistent with the observations include; a high 
surface temperature and low interior temperature during the very early lunar 
history; high near-surface radioactivity and relatively low radioactivity In 
the interior; a molten zone formed at or near the surface which gradually 
migrates downward with time. 

Lunar magnetic anomaly calculations show that large anomalies measured at 
some of the landing sites and above some points on the surface cannot be 
caused by mere basalts but are consistent with valley fillings of Cayley-like 
material with a remenant magnetization of about 2 x lO"^ emu/gm. 


We also speculate that the source of the magnetic field which must have been 
present in the early stages of lunar evolution could have been caused by a 
layer of conducting fluid at a depth of several hundred k i 1 acted 

to decouple the solid core from a solid crust. 
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This report was prepared by Earth Sciences Research Inc. and covers work car- 
ried out under contract NAS-9-11968 during the period June 1971 through August 
1973- The work was administered by the National Aeronautics Space Administration, 
Lyndon B, Johnson Space Center, under the technical direction of Dr. Paul W. 

Cast until his untimely death in 1973* Throughout the program close liason was 
maintained with Dr. Cast, Dr. David W. Strangway, Dr. Jeff Warner and 
Dr, Wolfe Gose. 

This report presents a summary of the findings of the study, much of which has 
already been incorporated into publications listed in the back, !t also in- 
cludes some unpublished material related to the thermal history of the moon. 
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INTRODUCTION 


i 


Prior to the Apollo landings much effort was devoted to predicting the possible 
thermal state of the lunar Interior. The investigations included development 
of computer programs to simulate the thermal development of the moon and thus 
provide a guide for the intelligent and systematic exploration of the lunar 
surface. Following the landings it was possible to eliminate many lunar his- 
tory models and to concentrate on refinement of those models which were found 
to be generally in accord with the evidence, obtained from the returned samples 
and other observations. The present authors had participated in one of the 
prelanding studies involving development of lunar thermal history simulation 
programs. Objectives of the research reported here were to employ actual lunar 
sample data in updating and modifying the existing computer programs in order 
to provide refined models of geological history of the moon and similar small 
planets and, subsequently to calculate thermal histories for specific parameters 
of Interest in establishing the evolution of such planetary bodies. 

The bulk of the effort was directed towards developing a computer program with 
the ability to handle such factors as: non-uniform initial composition; ar- 

bitrary initial temperature gradients; arbitrary initial distribution of radio- 
activity; effective latent heat of melting; effects of variable degrees of 
partial melting; effects of redistribution of heat sources; and effects of 
variable thermal conductivity, in addition, there were several minor modifications 
required such as modernizing the computer language, resolving some minor dis- 
crepancies in the programs and attempting to eliminate some of the known less 
realistic aspects of the computer model. 
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Following completion of the program mod i f icat ions , and as a result of the Apollo 
discoveries of a substantial lunar magnetic field, the study was extended to 

•s. 

Include model calculations on various aspects of the lunar magnetic field in- 
eluding anomalies to be expected above circular craters and anomalies in the 
neighborhood of several Apollo landing sites, 
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STUDIES OF THE THERMAL HISTORY OF THE MOON 


Early calculations by Urey, MacDonald and others were based primarily on the 
assumption of an initially uniform model of either ’’chondri t ic" or “terres- 
trial*' composition and assumed a thermal conductivity in which radiative 
transfer played a large role. As it was recognized that melting of the lunar 
interior could significantly Influence the thermal history, several attempts 
were made to take this into account, Fricker, et.al. (19^7) allowed for 
upward diffusion of radioactive elements and assumed that once complete melt- 
ing had taken place, convective transfer would dominate, McConnell, et.al, 
( 1967 ), on the other hand neglected convective transfer in the conventional 
sense, but allowed for transfer to the surface of radioactive elements with 
the melt, thus simulating diffusion. 


More recently new information has been obtained which demonstrates that better 
models are required. In the first place, analysis of lunar samples shows that 
an Initial uniform cbeml cal . compos i t i on is no longer tenable, and that neither 
the “chondritic** nor the “terrestrial" models is suitable for an' initial starting 
composition (Cast, 1971). Furthermore, recent experimental and theoretical 
studies show that radiative conductivity is unlikely to dominate heat transfer 
from the lunar interior (Aronson, et.al. 1967, 1970; Fukao, et.al. 1968). 
Accordingly, Wood (1970), Hays (1971) and Toksoz, et.al. (1971) have incorporated 
more real istic coniposi tions and thermal conductivities into their models and 
have attempted to allow for the redistribution of radioactivity at the time of 
me 1 1 i ng. 
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Even with these inprovements however there are no published calculations which 
intimately associate movement of the radioactive elements with movement of a 
melt, whose migration is controlled by reasonable mechanical and thermal 
considerations. Such a modei should be capable of predicting a wide variety of 
properties which could be compared with field observations. Some of these 
observations are discussed in more detail in the following section. 

Observational Constraints on Successful Models 

With the rapidly increasing knowledge of the geological, geophysical, and 
geochemical properties of the moon, the constraints imposed on successful models 
become more and more stringent. The most Important are those related to 
present interior temperatures, present heat flow composition, density and velocity 
distributions, strength of the interior, and thermal properties such as specific 
heat,- thermal conductivity, and heat of fusion. In addition, the model must 
predict times, rates, and volumes of extrusives consistent with those inferred 
from other studies. 


One of the most difficult factors to determine Is the present interior temperature. 

♦ • 

From magnetic observations, Sonnett, et.al. (1970) have inferred temperatures as 
low as 800® K at 500 kilometers and 1200^ K at 1000 kilometers below the surface. 
Such temperatures are most easily explained by very low uranium concentrations 
in the deep interior, as pointed out by Hays (1971). Tozer (1971) believes that 
the equilibrium temperature for convection by means of solid state creep would 
not be much different than that postulated by Sonnett. Since convection requires 
an increasing temperature as a function of depth, and the geochemical evidence 
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suggests that the moon was Initially much cooler at depth than near the 
surface, Tozer's convective transport mechanism is only likely to become 
effective at those times and in those regions where temperature Increases 
faster than the adiabatic gradiant. (n the calculations discussed below, the 
effects of solid convection are neglected, although we do not believe they should 
be dismissed. 

The presence of mascons provides additional evidence for low interior temperatures. 
Their maintenance would require enough strength down to 150 kilometers (Kaula, 

1969 ) to maintain stresses on the order of 100 bars over periods of more than 
io3 years. This Implies temperatures in the range SOO^llOO® K or less at these 
depths since the time that the lunar mare basalts were extruded 3.5 billion 
years ago. 

The only estimate of the lunar heat flow available at the time of the calculation 
(Langseth, 1970 was 3.3 microwatts per square centimeter. When considered in the 
light of magnetic and gravitational observations this value was unexpectably large, 
yet it is consistent with intense volcanic activity In the past and there is no 
reason at the present time to suggest that it is significantly in error. The pos- 
sibility exists that conditions at the Apollo 15 landing site are anomalous, but 
this can only be determined after several different sites have been investigated. 

Thus we have two apparently conflicting sets of constraints. The high heat flow 
and intense volcanic activity suggest extensive melting of the Interior both 
now and In the past. The gravl tat ional and magnetic interpretations seem to 
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require a moon which is currently relatively cool to depths of at least 200, 
and possibly 500 kilometers. Whether any model can s Imul taneous ly explain the 
conflicting observations remains to be seen. 

As our model allows for redistribution of radioactive elements, and provides an 
estimate of the average composition of the different portions of the moon as 
a function of time, the observed surface composition can then be compared with 
that predicted in order to establish the validity of any model. Unfortunately, 
current estimates of the overall surface composition are based on samples derived 
almost exclusively from the maria which is representative of less than half of 
the total lunar surface. However, the composition of the basalts generated in 
the interior must be consistent with the average basalts seen on the surface. 

This means that the uranium concentrat ion for the mare basalts generated about 
3.5 X 10^ years ago must be approximately 0.25 Ppm uranium and 500 ppm potassium. 

The. model should also predict an initial period of melting and extrusion whose 
age is greater than 4.4 x 10^ years and a second period of extrusion centered 
around 3.5 x 10^ years ago, which lasts for approximately 10^ years. Minor 
isolated volcanic activity throughout the reamlnder of the moon's history is 
permitted. 

The redistribution of major components implies a change in density distribution 
with depth which also must be reflected in the moment of inertia,} This requires 
a knowledge of specific volume of material as a function of pressure and 
temperature. At the present time we are assuming that the density of the low 
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melting phase both in its solid and liquid states is grams/cc and the 
high melting phase is 3.^ grams/cc. Because of this obvious overs impl i f ication 
of the density of the major phases, no attempt: has been made to calculate the 
moment of inertia. 

Geometrical, Compositional and Thermal Parameters 

As a first step in developing a more geologically complete model we modified 
our previous computer programs to bring them to a level of sophistication 
sufficient to incorporate the following features: 

1. A nonuniform initial composition in terms of fraction 
of low melting to high melting phase present and for 
variation in the uranium, potassium, and thorium contents 
as a function of depth. 

2. Partitioning of the radioactive elements between the melt 
and the solid phases, 

3. A cutoff value of melt which must be exceeded before 
magma can move to the surface. 

The model is in many ways similar to that described in a previous paper 
(McConnell, et.al. 1967)* As in the previous model, the moon is broken up 
into a n’umber of zones. For the calculations described below each zone was 
taken to be. 20 kilometers thick. For simplicity In comparing results with 
different zone thicknesses the radius was assumed to be 1700 kilometers. 

Three major compos i t ional phases whose concentrat ions can be specified as a 
function of depth are recognized: a high melting or “dunite" phasej a low 
melting or "basalt*' phase; and a magma phase. Three minor elements are also 
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recognized: potassium, uranium, and thorium. The potassium and uranium 
concentrations are specified independently as a function of depth. The uranium 
to thorium ratio is held constant at 3.7« 

The initial temperature, melting temperature and heat of fusion are specified 
as a function of depth, and the thermal conductivity is specified as a 
function of the local temperature. 

For most minerals the specific heat increases as a function of temperature, 
however, we have taken a mean value of 1.1 joules/gram *K as representative 
over the range of interest. The heat of. fusion also should be expected to vary 
with changes in the mineral assemblage present, again an average value, in this 
case 500 joules/gram has been chosen for the melting of the basaltic phase. 

The melting temperature has been represented by the solidus for the Apollo 
basalts as determined by Ringwood and Essene (1970). 

The two most uncertain thermal parameters In our calculations are the temperature 
distribution at the time of formation and the thermal conduct i vi ty . There 
js strong evidence that the surface of the moon was at or near melting at 
the time of formation and that the interior was relatively cooler, which has 
been attributed by some (Mitzutani, et.al., 1971; Wood, 1970 to differences In 
energy of accretion ^ a simple model for the Initial temperature distribution 
was taken where the temperature in '’K at radius r is approximated by the 
expression 

T (r) = fE/c +273 
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whe re ; 


E is the accretional energy in jouies/gram 
c is the specific heat 

f is the fraction of accretional energy retained. 

For the purpose of this calculation E has been taken to be the. energy of a 
particle arriving at the escape velocity appropriate to the particular radius, 

C to be 1*1 joules/gram, and f a constant .35, independent of depth. 

The relevant initial and melting temperatures are shown in Figure 1. 

Magma Migration 

The main difference between these calculations and those of other calculations 
Is the manner in which the movement of radioactive elements Is tied directly 
to the movement of melt. Under the assumption that infinitesimal amounts of 
melt will not move to the surface, the program accepts as input the number of 
grams of melt per cubic centimeter that must remain at a given depth after 
extrusion has taken place. The program simulates extrusive activity by removing 
any excess melt from a partially molten zone and adding It to the surface layer. 
Solid material, from above is then allowed to fall down and fill the space and 
the composition Is adjusted accordingly. 

For some of our early calculations liquid was allowed to remain In the zone 
above to block downward diffusion of radioactive elements concentrated in the 
melt. Unfortunately, this resulted in large concentrations of these elements 
at the bottom of the crust. After a while the concentration of heat sources 
there was so great that a continuous cycle of melting, extrusion, and remelting 
of the uppermost basaltic layers occurred. The resulting thermal profile 
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developed a large hump much like that appearing in Wood's (1970) calculations. 
Since this concentrat ion effect seemed unreal! st ic^ the procedure was altered 
so that solid was moved downward but the melt that would have been associated 
with it was tapped off and allowed to move to the surface. 

The effects described above occurred when the radioactive elements were assumed 
to concentrate entirely in the melt. The computer program also allows one to 
chose two other ways to distribute the radioactivity; either evenly betv;een the 
melt and solid basalt phases, or equally among the melt, solid basalt, and 
dunite phases. The results of theninal history calculations using the latter two 
options were described in considerable detail in McConnell, et.al. (19^7), but 
they only have been utilized during the current study for test purposes.. A 
flow chart for the entire thermal history program is shown in Figure 2. 

Results of Calculations 

A number of different models were tested during the course of^program development. 
None of these have yet been found to produce results in accord with all the 
observations. Because there is no simple linear relationship between the 
observations and the Input variables, the problem of finding a suitable model 
is still largely one of trail and error. Vie discuss below the results of 
several attempts to determine whether reasonable conditions, composition, and 
thermal properties can be expected to give rise to. two separate periods of 
volcanism; the first occurring between k.(> x 10^ and 4.1} x 10^ years B.P., and 
the second approxi mately 3*5 billion years B.P, 
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In all the models described below, the Initial composition was taken to be 
1,6 g/cc of the low melting or '*basalt“ phase and 1.7 g/cc of the high melting 
or "dunite" phase. The fraction of melt which' must be present before extrusion 
can take place was set at .8 grams/cc, approximately one quarter of the total 
volume. 

All models had a uranium concentration of .1 ppm at the surface which decreased 
to .018 ppm in the deep interior and a near surface potassium concentration of 
150 ppm, which increased to 350 ppm below. 600 km. (Figures 3» 5) GM-AA 

and are completely identical except that the former utilized the early 

version of the program discussed above which allowed a high concentration of 
radioactive melt to accumulate near the bottom of the crust, while GM-AB and 
GH“6 used the, later version which eliminated this effect. 

For GM-kA and GM-AB the thermal conductivity was taken to be .008 watts/cm “K. 
This value is very low for most basic and ultrabasic silicate rocks and roughly 
corresponds to the mean value of the synthetic lunar basalt reported by 
Murase G McBirney (1970)* It was chosen as a lower bound for reasonable 
conductivities in order to test whether melting could occur under the most 
favorable conditions in an initially unmelted moon. In Model GM-6 the 
conductivity was doubled to .015 watts/cm “K, a more typical value for basalts. 
In both cases changes of temperature with conductivity were neglected. 

Both versions of GM~^ predict volcanic activity beginning about x 10^ 
years B.P. and gradually decreasing until the present time. (Figure 6) 

The most recent activity Is considerably more intense for than for 

as a result of the continual overturn of the crustal layers caused by the 
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concentration of radioactivity at the base of the solid crust. The improvement 
in GM-4B resulting from abandoning this feature is obvious. Unfortunately, 
neither of the above models shows a distinct break between the early melting 
period and the later one.’ One possible explanation for this is the relatively 
smooth variation in radioactive element concentrations as a function of depth 
which results in a continuous downward migration of the melting zone. To 
test this hypothesis a much sharper gradation in the uranium concentration was 
utilized in Model GM-6 (Figure 5). We also Increased the thermal conductivity 
for this mode) to 0,015 watts/cm ®K, since other tests had shown that the 
extremely lov/ values were not necessary to induce early melting. As shown in 
Figure 6, GM-6 exhibits In a somewhat cleaner separation of the various periods 
of volcanic activity, but still not as distinct as indicated by field evidence. 
It Is not yet possible to say whether a better selection of parameters can 
improve the predictions or whether it will ultimately be necessary. to invoke 
other factors such as retention of the melt at depth until release is triggered 
by mare basin formation. 


Among the more successful predictions of the three models is the minor element 
composition of the mare basalts. The K, U, and K/U values for the upper 20 
kilometers at 3*5 x 10^ years B.P, are shown in Table 1, 


Table 1 


Model 

ppmK 

ppmU 

K/U 

GM-U 

A57 

.256 

1780 

GM-4B 

595 

.333 

1780 

GM-6 

38) 

.163 

23AO 
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AU are sufficiently close to the observed values to be considered reasonable. 


The present heat flow for the three models is 'shown In Table II. 


Table II 


Model 

Conduction 

Volcanic 

Total 

gm~4a 

1.03 

1.00 

2.03 


1.87 

.37 

2.24 

gm~6 

1.77 

.A6 

2.23 


As all show continuing volcanic activity at the present time the heat from 
the interior is thus transported via two mechanisms conduction and convection. 
The contribution of the volcanic heat to the flow measured by a lunar surface 
experiment is a strong function of the way the melt behaves when- it reaches 
the surface and is not predictable from the model. The total heat flow however 
is relatively insensitive to the transport mechanism and it is this which should 
probably be used to compare predictions and observations. We note that the 
predicted values for all three models are about thirty percent lower than 
that measured by the Apollo 15 heat flow experiment, suggesting that higher 
radioacti vi ty would provide a better fit. An alternative explanation is that 
a very high effective conductivity, perhaps due to solid state convection is 
causing the interior to cool rapidly at the present time. There is some doubt 
however as to v/hether this could provide a suitable explanation as the total 
heat flow is almost exactly equal to the present rate of production. 
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We turn now to the present temperature dfs tri but Ions (Figure 7). All models 
show Internal temperatures which are much higher than those predicted by 
Sonnett, et*al. While the near surface values could be reduced by increasing 
the conductivity and/or the extent to which radioactive elements are 
concentrated in the crust, we find It difficult to see any mechanl sm for 
lowering the temperatures In the deep Interior without a further drastic 
reduction in the original uranium concentration at depth. 
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Ill, LUNAR MAGNETIC ANOMALY CALCULATIONS 


During the Appolo program both surface magnetometers and subsatellite magneto- 
meters indicated the presence of magnetic field anomalies on the moon. In 
order to test various hypotheses regarding the origin of these anomalies a 
variety of calculations were carries out using the numerical techniques 
developed by Talv;anl ( 1365 )• 

We first tested the hypothesis that anomalies associated with lunar craters 
could be explained by uniformaly magnetized sheets of basalt, with holes, 
punched through them by impacts. Calculations were made of the horizontal 
and vertical field anomalies, as seen from a satellite passing directly over 
a crater at a height ^of 110 kilometers, for craters with radii of 50, 100, 
and 200 kilometers and thicknesses of 20, AO, and 100 kilometers. Three 
directions of magnetization were examined: vertical; horizontal parallel to 

the satellite path-; and horizontal perpendicular to the satellite path. The 
results of these calculations are shown in Appendix 1, 

To examine the surface anomalies at the Appolo 15 site a simplified map of 
the regions topography was drawn up and digitized (Figure 8). Because of the 
irregularity of the topography and the resulting problem of handling dis- 
connected contours, considerable time was devoted to developing appropriate 

techniques and running test programs. Calculations were then made using 

-6 

a magnetization of 2 X 10 emu/gm. which indicate an anomaly of approximately 
0.2 Y in the neighborhood of the LM. As a result of these calculations it 
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FIGURE 8 TOPOGRAPHIC, MODEL FOR MAGNETIC CALCULATIONS - HADLEY REGION 



v/as concluded that the anomaly at the landing site cannot be explained by 
rocks with magnetizations comparable to those sampled at the edge of the 
rille but can be explained by magnetizations which fall within the range 
of mare gasalt samples. This is discussed in more detail by Strangway et 
al. (1973). 

The largest surface anomalies measured during the Appolo program occur at 

the Appolo l6 site where 5 readings fell between 121 and 313 gammas. To test 

the hypothesis that these anomalies are associated with the Cayley formation 

breccia flow, which fills the valley and which is known to have strong remanant 

magnet ization, we computed the anomalies for a simplified model representing 

the edge of a basin filled with Cayley-like material. The results which 

-4 

indicate that a permanent magnet i zat i on of the order of 2 X 10 emu/gm. 
would account for the level of observed anomalies are reported in more detail 
in Strangway et al. (1973). 
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IV, SPECULATIONS ON THE RELATIONSHIPS BETWEEN THE THERMAL AND MAGNETIC 
HISTORY 


The apparent paradox regarding the thermal evolution of the moon at the 
time the heat flow calculations were carried out can be summarized as -■ 
-foUov/s: 

1. Paleomagnet i c studies of lunar basalts demonstrate that all 
such rocks returned to date were extruded at a time when the 
moon had a magnetic field. 

2. In order to have a magnetic field the moon is assumed to have 

had a molten metallic core. ‘ 

3. If the moon once had a deep molten core then It would still have 
one because of the low thermal conductivity. 

A. Present measurements of the moon's internal temperature based 

on solar wind and seismic observations indicate that the Interior 
is now re^l at I vely cool. 

5. Therefore the moon could never have had a molten core. 

The difficulty is to find a model warm enough to produce a mobile conducting 
layer and to melt basaltic rock and yet one which is cool enough in the deep 
interior to remain solid throughout the lunar history. 

Except for the magnetic field paradox most other observational constraints 
can be fitted by a model which involves an initially differentiated moon 
enriched in refractory elements near the surface. The geochemical evidence 
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for such a model has been summarized by Cast (1972). If such a differentiated 
moon had the nov%f widely accepted initial temperature distribution which was 
high near the surface and low in the Interior it most likely began melting 
near the surface and the molten zone migrated downward as time wnet on 
(McConnell and Cast 1972). Typical rates of downward migration are .3 km 
per million years, i.e. 10“? meters per second. 

The molten zone would tend to accumulate high density components including 
such low viscosity conducting melts as metallic sulfides and metals. The 
importance of such a low viscosity layer Is evident when one remembers that 
a molten core serves two purposes: first by being a highly conducting fluid 
it provides the physical environment necessary for generation of a magnetic 
field; secondly if the temperature gradient is superadiabati c convections 
will give rise to motions necessary for field generation. 

It is important to note that convection is just one way of developing fluid 
motions necessary for field generation. Presumably any other mechanism 
which can produce appropriate velocity distributions would be equally good. 

We show below that if during the period of orbital evolution the moon was 
divided into two regions, a solid outer shell and an inner core, separated 
by a layer of molten conducting liquid of the order of tnes of kilometers 
thick then the core would tend to rotate differentially with respect to the 
shell. If the magnetic Reynolds number for the flow in the decoupling zone 
were high enough then magnetic field generation could be expected even though 
the bulk of tfie core rema cool . 
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To demonstrate the feasibility of this mechanism let us take an oversimplified 
lunar mode) consisting of a sphere with a uniform interior core of radius r_ 
and density and an outer shell with sufficient gravitational dipole 

moment that remains oriented with one side aligned facing the earth. If we 
define (FI gure 9] ) 

^5 = angular velocity of shell relative to Inertial space 
“ angular velocity of core relative to inertial space 
W ' c ^ s 

then 



v/he re 

is the torque exerted on the core by the shell 
moment of inertia of the core 


For a molten layer of thickness h and viscosity r) the torque exerted 
on the core due to .the differential rotation o) can easily be shown to be 
approximately 





(jj 


( 2 ) 




( 3 ) 


Let us neglect the effects of changes in core size with time as a result of 
downward migration of the molten zone and concentrate on the relationship 
between the orbital changes and w . 
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Substituting the moment of inertia for the core 


Ic = 8 TT p ’ rc 


5/15 


and ( 2 ) into (1) yields the equation of motion 


w + 1 0 n 


i I ^ = 0 

h ) fij SJ^ 


which can be ' i nteg rated numerically if is knovm 


(A) 


( 5 ) 


For the sake of the present discussion let us assume that for the period 
with which we are concerned the moon v/as in an approximately circular orbit 
with the shell aligned towards the earth so that 


U 


s 


= 2.0 X 10^ 


where 


( 6 ) 


G is the gravi tat ional constant = 6.67 X 10 ^^Nm/kg^ 

24 

M is the mass of the earth = 5.98 X 10 kg 
R is the radius of the mQon ’ s orbit 


If the torque acting on the moon Is related to the radius by an expression 
of the form 


L = k R”" = ^ (m n) 
dt 


v/he re 

k is some constant 
m is the mass of the moon 

is the orbital angular velocity 

33 . 
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It can easily be shown that 


= -3 R = (8) 

fij 2 R ltS+r/2) lt2-t'j) 

Here ti and t 2 are the times at which the moon was at radii Rj and 
R 2 respectively. 

Choosing 

a = 6 appropriate for torques due to earth-moon tidal interaction 
R| = 20 earth rad i i 
^2 “ 60 earth rad t i 
tz“ti = ^.5 X 10^ years 
Pe = 3.3 X 10^ kg/m^ 

n = 5 X 10 ^ kg/m sec, typical of molten Iron 

o g 

w = (0 = 0 at t ^ -^.5 X 10 years 

we then integrate (5) to determine oi(t) for various layer thicknesses h. 
Having obtained w we are now in a position to determine when, if at all, 
conditions were favorable for magnetic field generation. 

One normally assumes that a combination of turbulent flow and a magnetic 
Reynolds number much greater than 1 is required for dynamo action. 

For large Reynolds numbers flow between rotating cylinders is unstable when 

'^s(rc + h)^ < S5^(r^.)^ ( 9 ) 
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3k 


i^arth 


ctehces 



(Landau and Lifshitz 1959 p.l08). If we assume that the same rule Is 
valid for flow between spheres, substitute Qc " fls ~ w (8) • and 

rearrange v;e get as a criterion for instability 


0 = w 25. + > 1 

^ --C ^ 

The magnetic Reynolds number is defined as 
= ^vpcr 

where 


5..' is a character i St i c length 

V is a characteristic velocity 

V is the magnetic permeability 

o is the electrical conductivity 


(to) 


(n) 


For this particular problem we may take 

|i = 4ir X 10" webers/m^ 
o = 3 X 10^ mhos/m and 
i^v = r^ - r^) 

whi ch s i mpl i f I es to 

Re = h^ - whr_) (12) 

jt; s c 

For the moon dynamo the ratio of viscous Reynolds number to the magnetic 
Reynolds number: P/hUO >> 1 
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Computed values of B as a function of time for h = 10 km and h = 20 km 
are shovvn if Figure 10. Values of Re^ are of the order of 10^ during the 
unstable period. Although h = 20 km would appear to give the best fit to 
the observations, it should also be pointed out that lower viscosities, 
which might be more appropriate for sul fi des, would have the same effect. 

One must also remember that (5) is valid only for laminar flow and that the 
onset of instability will give rise to increased viscous drag and drag 
induced by the build up of the magnetic fjeld. The net result of this 
will probably be to stabilize the system at a value of 3 approximately 
equal to | rather than allowing it to overshoot as shown .in Figure iQ., 

V/e therefore conclude that if a conducting layer with low enough viscosity 
to effectively decouple the core and outer shell is formed beneath the 
surface of the moon then conditions favorable to magnetic field generation 
will arise as a result of the normal orbital evolution and that much of the 
reason for the hot versus cold moon controversy no longer exists. 
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V.- SUMMARY AND CONCLUSIONS 


On the basis of the calculations discussed in the preceeding chapters we conclude 
that a realistic computer simulation of the evolution of the moon can be used 
to predict many more observables than heat flow and internal composition. These 
observables include major and minor element concentrations, thickness of the 
lunar "crust", intensity of volcanic activity as a function of time, etc. r 

The models which are most consistent with the observations include the follow- 
ing features: 

1. A high surface temperature and low interior temperature 
during the very early lunar history; 

2. High near surface radioactivity and relatively low radio- 
activity in the interior; 

3. A molten zone formed at or near the surface which gradually 
migrates downward with time. 

The lunar magnetic anomaly calculations demonstrate that the large anomal ies 
measured at some of the landing sites and above some points on the surface 
cannot be caused by mare basalts but are consistant with valley fillings of 
Cayley-like material with a remnant magnetization of about 2 X 10 ^ emu/gm. 

We also speculate that the source of the magnetic field which must have been 
present in the early stages of lunar evolution could have been caused by a layer 
of conducting fluid at a depth of several hundred kilometers which acted to de- 
couple the solid core from a solid crust, 
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MAGNETIC ANOMALIES OVER MODEL LUNAR CRATERS 


It has been postulated that magnetic anomalies associated with lunar craters 
result from the edge effects from around a hole punched out of a uniformly 
magnetized basaltic crust. If this is the case the measured anomalies can 
be used to estimate the thickness of the crust and the direction of 
magnetization at the time of implacement. These values in turn impose important 
constraints on the thermal and volcanic history of the moon. 

The Mode 1 s 


The anomaly due to uniformly magnetized sheet with a hole punched in it is 
assumed to be equivalent to that of a disc of the same size as the hole with 
equal but opposite magnetization. For the purpose of the calculations such a 
disc is approximated by a regular sixteen sided polygon and the anomaly 
calculated using the method of Taiwan! (1965). 

We have calculated anomalies for a variety of crater sizes, slab thicknesses, 
and magnetization directions in order to determine the main features of such 
anomaly patterns. For a coordinate system with origin at the top of the 
slab in the center of the crater the anomaly field components AX, AY, and AZ 
were calculated along the positive x axis out to the distance of AGO kilometers 
for the following models: 


Satellite hei ght 
Crater radius 
Slab thickness 
Magnet i zat i jon 


110 ki 1 ometers 
50, 100, and 200 kilometers 
20, AO, and 100 kilometers 
No susceptibility effects 
Remnant Intensity 10“^ emu/cc 

Remnant direction - horizontal parallel to satellite path 

~ horizontal perpendicular to satellite path 
- vertical 


Presentation of Results 


The results of the calculations are given in the accompanying tables and figures. 

For each crater three tables have been calculated, one for each direction of 
magnet i zat i on. With suitable transformations these can be used in conjunction 
with a desk electronic calculator to estimate the anomaly for any other field 
orientation or position over the crater. 

Two graphs are also shown. The first presents the anomaly for the two 
horizontal magnetization directions as follows: 

1. For the slab with the magnet i zat i on J parallel to the x axis 

the aX and aZ components are indicated by — and f 
respectively. aY = 0 and is not shown. Y 

2. When the magnetization Is parallel to the y axis only the aY 
component is non-zero. This is indicated by Q . 

/ 



The second graph shows the anomalies for vertical magnetization with: 
AX indicated by ; 

AY wi)ich is equal to zero, Indicated by Q ; 

AZ indicated by | 
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'-0’28l” 

0.000 

o.ooc 

0,000 

80.00 

0,0 

0.000 

-0.203 

oiboo 

100 . no 

0.0 

”0. QOO 

- 0 * 1 n 

-n . non 

120.00 

0*. 0 

o.noo 

-0. 132 

0.000 

j,40. on 

0.0 

- 0 . non 

“O.ins 

- 0 . non 

IGO. 00 

0.0 

“0.000 

-0.0^4 

“0.000 

leo.oo 

0^ 

-J).noo 

-0.008 

_D.i)00 . 

200.00 

0.0 

-o.noo 

-0.035 

0.000 

220.00 

o.n 

0.000 

-0.045 

0.000 

240.00. 

0.0 

-o.noo 

-0.037 

0 . non 

2G0. 00 

0.0 

“ 0 . non 

-0.031 

-0. non 

280.00 

0.0 

0.000 

-0.02C 

0.000 

300.00 

0.0 

-o.noo 

-0.022 

0.000 

320.00 

0 , 0.. 

0.000 

-0.019 

0.000 

340.00 ■ 

0.0 . 

b.ooo 

-O.OIB 

b.ooo 

560.00 

b. b 

0.000 

'"-0.014 

n . nob ’ 

380. 00 

0.0 

-0.000 

. ■ -0.012 

0.000 

400. 00 

0.0 

0.000 

-C. 010 

„_-D. JIOJI, 






SctSiices 




SATELLITE HT. = 110 Km 

RADIUS = 50 Km 

DEPTH = 100 Km 

MAGNETIZATION 

INTENSITY = I0“® emu 
DECLINATION 0.0 deg 

INCLINATION = 90^0 deg 



SATELLITE 

POSITION 


ANOMALY 


K 


AX 

AY 

AZ 

0.0 

. 0.0 

- o . ooo ' 

0 . 0 0 0 

0 . 633 - 

20.00 

0 . 0 

~(). 152 

0 . 000 

0,3 9 b 

40.00 

0.0 

tD .220 

0.000 

O .' i 0 (> 

GO . 00 

0 . 0 

- 0'.'2 74 

0^,0 ao 

0 ._5 0 8 

80.00 

0.0 ' 

- i ).271 

0.000 

0 . 24 G 

100. no 

0.0 

< - 0.239 

0.000 

0.150 

120.00 

0.0 

- 0.197 

0.000 

0.025 

. 140.00 

0.0 : 

-( LJ .57 

- 0 . 000 

0.043 

IGO . 00 

0.0 

- 0.123 

- 0.000 

0.019 

180.00 

0.0 

- 0 . 096 

0 . non 

o ._ on 5 

200.00 

0.0 

- 0.075 

0.000 

- 0.003 

220.00 

0.0 

- O .059 

0.000 

- 0.007 

240. 00 

0.0 

- 0.047 

0.000 

- 0.009 

200.00 

0.0 

- 0.037 

- 0. 000 

- o.oin 

280.00 

0.0 

- 0.030 

0.000 

- 0.010 

300.00 

0 . n . 

- 0.025 

0.000 

- 0.010 

320.00 

0.0 

- 0.020 

0.000 

- 0.009 

540.00 

0.0 ■ 

- 0.017 

0.000 

- 0.000 

300,00 

0.0 

- 0.014 

0.000 

- 0.008 

380.00 

0,0 

. 

n.oon 

- 0 . 007 

400.00 

0.0 

- 0.010 

- 0.000 

— --La-.V.* «|— fxr -UM . 

- 0.007 - 




X 


SATELLITE HT. 

RADIUS 

DEPTH 


= no Km 
= I 00 Km 
= 20 Km 


MAGNETIZATION 

INTENSITY = I0"° emu 

DECLINATION >= 0.0 deg 

INCLINATION = 0.0 deg 


SATELLITE POSITION 



ANOMALY 


% 


X 

AX 

AY 

£ 

,z 

• 0. 

, 0 

0.0 ; 

-0.167 

0 . onn 

_Q. 

QOO 

20. 

.00 

0.0 . 

-O.lGl 

0.000 

-0. 

040 

40. 

. 00 

0.0 

-0.143 

-0.000 

-0. 

005 

■ 60. 

.00 

0.0 

-0.113 

0.000 

-0. 

132 

P.O . 

. 00 

0 . 0 

-0.075 

o.onn 

-0. 

157 

lo'b. 

. 00 

0.0 

-0.035 

0.000 

-0. 

164 

120, 

,00 

0.0 

0.000 

n, 000 

-0. 

156 

140 , 

. 00 

d. 0 

0.026 

0.000 

-0. 

137 

160. 

, 00 

0.0 

0.041 

-0.000 . 

-0, 

114 

1 0 0 1 

, 00 
, 00 

0 . d 

0.048 

o.iJoo 

-0. 

091 

200, 

0.0 

0.048 

0.000 

-0. 

071 

220. 

, 00 

0.0 

0.046 

-0.000 

-0. 

056 

043 

240. 

.00 

0,0 

0.042 

0,000 

-0. 

26 0. 

, 00 

d . b 

0.037 

0. 000 

-0. 

034 

28 0-. 

_00_ _ 

0.0 

0.033 

0.000 

“0. 

027 

300 . 

00 

d. d 

0.029 

0.000 

-0, 

021 

320. 

00 

0. 0 

0.025 . 

0.000 

-0. 

017 

340. 

00 

0.0 < 

0,022 

0, 000 

-G* 

01 h 

3 b 0 . 

00 

Ol^O 

0.019 

0.000 

-0.011 

380 . 
400. 

00 
qq _ 

0 . d 
._,_0_^0 

0.017 
- 0.015 

0.000 

-0.000 

-0. 

009 ■ 
008. . • 


^arth . 


/^^csearth. 



X 


SATELLITE HT. 

RADIUS 

DEPTH 


110 Km 
rOOKm 
20 Km 


MAGNETIZATION 


INTENSITY 

DECLINATION 

INCLINATION 


= 10"® emu 
= 90.0 deg 
= 0.0 deg 



CRATER 


SATELLITE POSITION 


ANOMALY 


X 

X 

AX 

AY 

AZ 

0.000 
. 0.000 

0. 000 
n Cl n n 

0. 0 
20.00 

4o;‘oo' 

CD.A0 

RO.OO 
100. on 

o.d 

0.0 

* 0.0 

_...a,.jQ 

n. 0 
0.0 

0.0 00 
0.000 
-0.000 
n 0 00 

-0.167 

-0.165 

-0.159 

-0-1 ii n 

oiobo 

-0.000 

-0.156 

-0.120 

0.000 

0*000 

120.00 
140. on 

O'. 0 
0.0 

0.000 

0.000 

-0.102 

-0*086 

o.coo 
- 0 . non 

ICO. 00 

0.0 

-0.000 

-0.071 

o’. 000 

180.00 

0.0 

0.000 

-0.058 

“0. 000 

• 200.00 

6.0 

0.000 

-0.047 

-0. 000 

220.00 

0.0 

-0.000 

-0.059 

-0.000 

240.00 

0. 0 

0.000 

-0.032 

0. 000 

2G0.Q0 ■ 

0.0 ■ • 

0.000 

-0.026 

0.000 

280. 00 

0.0 

0.000 

“0, 022 

0.000 

300. 00 

0.0 

0.000 

-0.018 

-0. 000' 

32"0. bo 

b. 0 

0.000 

-0.016 

0.000 

340.00 

0.0 

0.000 

. -0.013 

o.coo 

560.00 

6.0 

0.000 

-0.011 

0.000 

380. 00 

0.0 

0. 000 

-0.010 

0 . non 

400. 00 

0.0 

-0.000 

-0.009 

0.000 


ts ar til 


S’cisncis 





SATELLITE HT. = liO Km 

RADIUS = lOOKm 

DEPTH = 20 Km 


MAGNETIZATION 

INTENSITY =10"® emu 

DECLINATION = 0.0 deg 

INCLINATION =90.0 deg 


SATELLITE POSITION 


X 



CRATER 

ANOMALY 


K 

X 

AX 

AY 

AZ 

■ 0, 0 ■ ■ 

0.0 

~o. non 

0.000 

0. 33 5 

20,00 

0.0 

' - 0 , n 9 

0 . non 

0.327 

fjO, 00 

0.0 

“0.095 

0,000 

0.302 

GO. 00 

n.o 

.“.a...L3 2. ... 

.0.000 

jD.2 63 

80.00 

0.0 

“0.157 

0.000 

0.211 

.too. 00 

0.0 

“0.164 

0. 000 

0.155 

120.00 

0.0 

“0.15G 

0 . ono 

0.102 

■ IkO, 00 

0.0 

“0.137 

0.000 

0. 060 

IGO. 00 

0.0 

“0.114 

0.000 

0.030 

180. 00 

0.0 

-0..0 9T 

Q. aoo 

0.010 

200. 00 

0.0 

“0.071 

“0. 000 

“0.001 

220. 00 

0.0 

“0.056 

-0 J300 

“0.00 7 

2^f0. 00 

0.0 

“0.045 

0.000 

“0.010 

2G0. 00 

0.0 

“0.034 

0. OOP 


280. 00 

0.0 

“0.027 

0. 000 

“0.011 • 

300.00 \ 

0.0 

“0.021 

..Q...00J} 

“0. 010 

‘320. 00 

0.0 

“0.017 

0.000 

“0.010 

340. 00 

0.0 

„„-0.014 

Q^OOO 

“0. 009 

360. 00 

0. b 

“0.011 

0.000 

“0.008 

380.00 

n.o 

“0. 009 

0 . non 

“0.007 

400.00 

0.0 

‘ “0.008 

0.000 

“0.007 


^arth Sc(Chc&% '^cuarch, 




f- . . . 



SATELLITE 

RADIUS 

DEPTH 


HT, = HO Km 
= 100 Km 
= 40 Km 


X 


MAGNETIZATION 

INTENSITY 
DECLINATION 
INCLINATION 


= 10"® emu 
<= 0.0 deg 

= 0.0 dea 



CRATER 


SATELLITE POSITION 


ANOMALY 


X 

X 

AX 

AY 

AZ 

0.0 

!L 0 

“ 0.309 

'o.ooo 

0.000 

/ O . 00 

0.0 

“ 0. 298 

0.000 

.“ 0.090 

40 . 00 

0,0 

“ 0.265 

“ 0.000 

- 0.175 

GO . 00 

0.0 

" 0^211 

' d . 000 

“ 0.242 

80.00 

0.0 . 

- 0.143 

0.000 

“ 0 . 2 < S 7 

100.00 

O.'O • 

“ 0.069 

0.000 

- 0.501 

120.00 

0.0 

- 0.004 

0.000 

“ 0.288 

140.00 

0.0 

0.045 

0.000 

“ 0.2 54 

’ 160.00 

0,0 

0.071 

- 0.000 

- 0.213 

180.00 

0.0 

0.084 

o.ono 

“ 0.172 

. 200. 00 ■ 

0.0 

0 . 087 ' 

0 . 00.0 

- o ,, a 3 j . 

220.00 

0.0 

0.083 

- 0.000 

“ 0. 108 

240.00 

0.0 

0.077 

- 0.000 

- 0.085 

2 G 0. 00 

0.0 

0.069 

0.000 

- 0.067 

280 . 00 

0.0 

0 . 062 

, 0.000 

- 0, 053 

300 . 00 

0.0 

0.054 

0 . 0.0 0 

“ 0.043 

320.00 . 

0.0 

0.04 8 

“ 0.000 

- 0.034 

5 40.0 0 

0.0 

0.04 2 

0.000 

“ 0.028 

3 G 0. 00 

0.0 

0.037 

0.000 

“ 0.025 

380.00 

0.0 

0.033 

0.000 

“ 0.019 • 

400.00 

0.0 

0.029 

“ 0. 000 

“ 0.016 


t:^arth 


G. 


cicnces 


es&arch, 



SATELLITE HT, 

RADIUS 

DEPTH 


= no Km 
= 100 Km 
= 40 Km 


MAGNETIZATION 


INTENSITY = 10"^ emu 

DECLINATION = 90.0 deg 

INCLINATION = 0.0 deg 


X 



CRATER 


SATELLITE POSITION ANOMALY 


X 


X 

AX 

AY 

AZ 

0 . 

O "" 

o : 0 

. 0.000 

-0.309' 

0.000 

20. 

00 

0 . 0 i 

i 0.000 

-0.305 

0.000 

"40 : 

00 

0 ; 0 

' -0.000 

-0. 204 

0.000 

60. 

00 

0.0 

: 0.000 

-0.276 

0. 000 . 

"80. 

on 

0,0 , 

0.000 

-0. 2,51 

0.000 

ino . 

00 

0.0 

: 0.000 

-0.222 

0.000 

120. 

00 

0.0 


r . 0 .... 131 _„„ 

. 0 ._ 0 - 0 ,Q ■ 

140 . 

00 

0.0 

0.000 

- 0.161 

- 0.000 

IGO ". 

OO " 

0 ;o ■ 

- 0.000 

- D ^._ L 33 ^ 

0.000 

1 H 0 . 

00 

0.0 

. 0.000 

- 0.110 

0.000 

200 . 

00 

0.0 

n . 000 " ' 

- 0. 090 

- C.noo 

220 . 

00 

0.0 . 

- 0.000 

- 0.074 

“ 0.000 

2 . 40 . 

0 0 

0.0 

0.000 

- 0, 061 

0.000 

2 G 0 , 

00 

0.0 

‘ 0 . 000 , 

- 0.051 

0.000 

23 0 . 

OG 

0.0 

;• 0.000 

- 0.045 . 

o.noo 

300 . 

00 

0 . 0 


- 0.0 36 " 

0.000 

320 . 

00 

0.0 

^ 0.000 

- 0.031 

0.000 

54 9 . 
3 G 0 . 

00 

o¥' 

0.0 

o.'o" 

0.000 

“ 0.026 

0.000 

3 30 . 

00 

0. 0 

OvO ^ jO 

0 -.- 02 - 3 — 

GrOOO . 

Too'; 

00 

'oTo 

0 . 000 
-- O ... Q . OJ ) 

-0.019 
r_C.JDlZ_ 

o.noo ■ 




Sciences 


i^c%&afch 



SATELLITE HT. 

RADIUS 

DEPTH 


= 110 Km 
= 100 Km 
= 40 Km 


MAGNETIZATION 

INTENSITY = I0“® emu 

DECLINATION = 0.0 deg 

INCLINATION = 90.0 deg 


X 



CRATER 


SATELLITE POSITION 


ANOMALY 


X 

_y_ 

AX 

AY 

AZ 

a. 0 

0.0 

^ 0*. 0,0,0. 


Q ...6_I8, 

20. 00 

0.0 

“0.090 

0.000 

0.60 3 

40, 00 

0. 0 , 

~.0...,,1.Z3.- 

_ „Q,...,0C0 

Q._i3.9,. 

60, 00 

0. 0 ' 

“0.242 

0. 000 

0 . /{ 8 7 

no. 00 

0. 0 



Q . 3 94 ‘ 

100. 00 

0,0 ' 

' “0.501 

0,000 

0.292 

UJ).J1SL 

0.. 0 

“0.,.,28 8,._.^ 

__..0..._OJ3,0 

0.196,; 

140.00 

0.0 

“0.254 " 

0.000 

0.118 

160. 00 

0.0 

“0.213 

0.000 

0,062' 

100, 00 

0.0 

“0.172 

0.000 

0.026 . 

200.00 

0.0 i 

-. ~A..U2 

. -:.a.ii.ao 

.. 0. OO’t 

220.00 

0.0 

“0 . 108 

-0.000 

“0.009 

^240. on 

0. 0 

“0.085 

0,000 

“_0,..015 

260, 00 

o.c 

“0.0G7 

0.000 

“0.018 

..280,00 

: 

..„~0, 053 

.0,.0„00 


300. 00 

0.0 

“0,043 

0.000 

“0,018 

320.00 

0,0 


Q_.iL0i) ' 

-0.017 

340.00 

0.0 

“0.028 

0.000 

“0.016 

.360. -00 

0,0 

. “0.023- 

0. 000 

“0.015 

380. 00 

0,0 

. “0. 019 . 

0.000 

“0.014 

40 0,. 00 



-r..Q.._ai,6. 

.a,„ao„o. 



tZiirtli Sciences J^csetrrch, "Vne, 



SATELLITE HT. 

RADIUS 

DEPTH 


X 


h 


= no Km 
= I 00 Km 
= 100 Km 


MAGNETIZATION 


INTENSITY 

DECLINATION 

INCLINATION 


= I0“® emu 
•= 0.0 deg 

= 0.0 deg 



50 Km 1 

■' ' — I 

CRATER 


SATELLITE 

POSITION 


ANOMALY 


X 

X 

AX 

AY 

AZ 

o \ 0 

0.0 

- 0.721 

0.000 

0.000 

20 . 00 

0.0 

- 0 . G 96 

0.000 

- 0.209 

on 

0.0 

- 0.620 

- 0.000 

“ 0.403 

^ GO . 00 

0.0 

- 0.496 

0.000 

- 0.564 

00 . 00 


- 0.335 

o.oon 

- 0.669 

100. 00 

0.0 ■ 

- 0.163 

0.000 

“ 0.704 

120. 00 

O '. 0 • 

- 0.011 

0 . ono 

- 0 . G 70 

•ino. 00 

0.0 

0.099 

0.000 

“ 0.590 

IGO.OO 

0.0 

0.163 

- 0.000 

- 0.493 

1 K 0 , 0 0 

0.0 ■ 

0.191 

0.000 

“ 0.399 

200. 00 

0.0 

0.195 

“ 0.000 

- 0.319 

22 0 . Of ) 

0.0 

0.187 

- 0.000 

“ 0.253 

240. 00 

0.0 

0.172 

- 0.000 

- 0.202 

260.00 

0.0 

0.156 

0.000 

- 0.161 

280. 00 

0.0 

0 . 140 

0.000 

- 0.130 

300. 00 

0.0 

0.124 

0.000 

- 0.106 

320.00 

0.0 

0.111 

0.000 

“ 0.086 

3 ' 40 '. 00 

0.0 

0.098 

0.000 

- 0.071 

360. 00 

0 . 0 ; 

0 . CS 7 

0.000 

- 0.059 

3 8 {) » 00 

0.0 

0 . 07 S 

0.000 

- 0.049 ' . 

400.00 

0.0 

• 0.069 

- 0. 000 

- 0.041 



Stisnces 






nc, 



SATELLITE HT. 

RADIUS 

DEPTH 


= no Km 
= 100 Km 
= 100 Km 


MAGNETIZATION 


INTENSITY 

DECLINATION 

INCLINATION 


= I0~® emu 
E 0.0 deg 
= 90.0 deg 


SATELLITE POSITION 


X 



CRATER 


ANOMALY 


X 

X 

AX 

AY 

AZ 

0.0 
20. 00 

0.0 
0. 0 

0.000 

-0.200 

0 . ono 
0 0 Q 

1 . '(2 
1 ^408 

40.00 

00.00 

0.0 i 
0 . JL 

-0.f|03 

-0..,5_64....... 

0.000 
. Q...a.o.a. 

1.307 

1.141 

80.00 
■ 100.00 
120.00 
140.00 

0.0 
0.0 
0.0 
0. n 

“0.669 

,,70jf „ 

“0.G70 

“0.590 

0.000 



0 . non 
o.noo 

0.924 



0.459 -■ 

0.278 

160. 00 
180. 00 

0.0 

“0.493 

0.000 
n . non 

0.150 
0^0 6 18 

200. O'O 
220. OO 

0.0 

0.0 

“0.319 

“0.253 

-n. 000 
“0.000 

0.018 
“0 . m 0 

240.00 
2G0. 00 

0.0 

0.0 

“0.202 

“0.161 

0.000 

0.000 

-0.026 
- 0 . 033 

230.00 

0.0 

-0.130 

0.000 

“0 . 036 

. 300. 00 

0.0 

“0. 106 

O.OOQ 

- 0 

320.00 

0. 0 

“0. 086 

c.ooo 

“0. 036 

3U0. 00 

0.0 

.-O.._0 7_l 

“0.000 

~0._Q.3fL 

300.00 

0.0 

-0.059 

0,000 

-0.052 

_3_B^00 

0. 0 

QJ1l9 , ■ 

-n. nno 

-0.029 

f}0~0. 00 

0.0 

-0.041 

0.000 

“0.027 





= 110 Km 
= 200 Km 
= 20 Km 

MAGNETIZATION 


SATELLITE HT. 

RADIUS 

DEPTH 


INTENSITY =10"® emu 

DECLINATION = 90.0 deg 

INCLINATION = 0.0 deg 


X 



CRATER 


SATELLITE POSITION 


ANOMALY 


X 


y 


AX AY A Z 


A A 

\) ^ \) 

20.00 

A A 

\i . l-> i 

■0.0 

40T0'0~^ 


CO. 00 

0.0 

CO ”.“00 

tlrO ■' 

100. 00 
— 1-A-A fVfl 

0.0 

A A 


IfjO.OO 0.0 

“1-COtO-O OtO— ■ 

120.00 0.0 


“2tr0T'0t) — 

220.00 

0 , 0 — 

0.0 ! 

.-^frOrdO 

2CO.OO 

0 . O""" 
0.0 ■ 

~2ot) ."OO' 
300. 00 

0 . 0 
0. 0 , 

• 320r0t) 
340.00 

0 ♦ 0 
0. 0 

“3"6urtj'u 

380.00 

0.0 . 


-OOOrOO OvP — 


nAA „A OAA „n nno 

'J 4 \J 0 U . i V./ V 

- 0.000 - 0.200 0.000 

'-OTtrotr- "-"0-.'200 o-rooo- 

0.000 -0.109 -0.000 

-DtOOD 

-n.ODO -0.194 0.000 

— 0; - i - e 9 “O^-e-0- 

0.000 -0.123 0.000 

0".-i-7-4 -&T-00O 

0.000 -0.1G2 0.000 

-OTi)^-0— — -G-.-HrS OttO 0 0 


0.000 

HhrO-5-e^ 

0.000 

-0-.ilOO- 

‘ 0.000 

- 0.-000 

O.COO 

- , /r\ f\ 

u ; V u u 


0 , 000 

-O-.-OdO- 


-0.132 

-o! 101 


-0.074 
“Gv-Olrir' 
-0.055 
•— GttOV?- 
-0.041 
— G-i-03-5- 


0.000 
-O^TdO'O- 
0.000 
-.-IHKI- 
0 . 000 
0.H3-OO- 
0.000 


-OvO-Od- 

0.000 


^arfh 3cie/]ccs H.Cr^ecirch,'^nc. 



SATELLITE HT. 

RADIUS 

DEPTH 


= no Km 
-200 Km 
= 20 Km 


MAGNETIZATION 


INTENSITY = 10"® emu 

DECLINATION = 0.0 deg 

INCLINATION = 90.0 deg 


X 



SATELLITE POSITION 


ANOMALY 


X 


y 


& X AY A Z 


; 

i 20.00 0.0 ; 

■'■“'0 0.00 

GO. 00 0.0 ■ 

^-so.00 OvO-"-". 

100.00 0.0 

"-1'207D0”" — ^OrO — 
140.00 0.0 

“IGOvOO OtO"— 

180.00 0.0 

“2 0 0; 00 0-."0 — 

220.00 0.0 

“240t00 0-rQ — 

2G0.00 0.0 

-•280v00 OrO — 

500.00 0.0 

--5 2 0.- 00 o-.-O— • 

540.00 0.0 

-360-00 0-0 — 

580.00 0.0 
~4 0 0.- 0 0-— 0-.- 0 — 


—O-T-nft-O -I 

! -O. 027 0 . 000 0.400 

—w07'0 n 5 0 rO 0 Q- ' 0':-5 0 8-- 

-0.084 0.000 0.595 

“0 71-16 00 Ot 5 8 4 

-0.150 0.000 0.3CO 

— "--0;-l-8 5" 0-rfH) 0""~- -e-r5fr6 

-0.218 0.000 0.512 

0 t- 2 4-6 OvO 0 0 : 0-2 G'5 • 

-0.261 0.000 0.208 

0-.--2-6 0 0-.-0 0 0 0-.“l 4-5 

-0.242 0.000 0.087 

e-.;“2-i-5 Gr-OilO- 0— O'^rl 

“0.178 0.000 0,003 

0^43 Q„„ 0 .. 3^„3 — _j 

-0.117 0.000 “0.024 

0 9 4___ Q 0 2-9— — - 

“0.075 O.OGO -0.051 

0%-G frO — CGr-O 0 ^ 

. -0.049 0.000 “0.029 

— ^--.O-.-0-5-9 0-.r-0.00 ^^0-,-0-2-Z ‘ 
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SATELLITE 

RADIUS 

DEPTH 


HT. = HO Km 
= 200 Km 
= 40 Km 


X 


MAGNETIZATION 

INTENSITY 
DECLINATION 
INCLINATION 


= 10*“^ emu 
= 0,0 deg 
= 0.0 deg 



SATELLITE POSITION 


ANOMALY 


AX 


AY 


AZ 


0,0 

„20. 00 

4 0.00 

-_G-0-.-0-0 

80.00 

100 -. 00 

120.00 

14 0,-00 — 


0.0 
-O-.-O— 
0.0 
- 0 ^ 0 - 
0.0 
0 .- 0 - 
0.0 
-OvO 


-30 OtOO- 

320.00 

540.00 
3GO.OO 
380.00- 

400.00 


-fh-0 — 

0.0 

■ 0.0 — 

0.0 

-O-.-O — . 

0.0 : 


-0.584 

-0.379 

-■*'0-'*-3-7-0- 


0.356 
■0-.--3 33- 
0.207 


IGO.OO 

_r» ft.A— A A 

0.0 

A — p 

-0.174 

J. O vj . \J \J 
200. 00 

0.0 

0,000 

- 2 2 0,-00 — 

— OvO 

-OvO-78- 

, 240.00 

O.D 

0.135 

-260 rOO--- 

— OvO 

— ^-0vl68- 

280.00 

0.0 

0.182 


0.173 

-0-.-160- 

0.146 

.O-vl-32- 

0.119 


0.000 
-O-.-OOa- 
0 . 000 
-O-r-O-O-O- 


0. 000 
-G.-000- 
0.000 
- 0 . ooe- 
0.000 

c.ooo 

0“i"0 O’O 
- 0.000 
0“0 o-o— 

0 . 000 

* ■■ Ml — .-I—— 

0.000 

0-000 

0.000 

0,-000 — 

0 . 000 


“ 0.000 
“0. lOG 
-0.223 

— Ov2-8G 

-0.351 

-Ov-4-1-3- 

“0.463 


-0.489 
—O-i-4 57- 
“0.403 

“0. 281 
-^-r-2-2-8- 

“0.184 
— 0il4G- 
- 0.121 

--..0 ..^0 9-8-.. 
-0.081 
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SATELLITE HT. 

RADIUS 

DEPTH 


X 


= 110 Km 
= 200 Km 
= 40 Km 

MAGNETIZATION 

INTENSITY = 10“^ emu 

DECLINATION = 90.0 deg 

INCLINATION = 0.0 deg 


SATELLITE POSITION 




Km 



CRATER 


ANOMALY 


X 


y 


AX AY A Z 


Q_0 GtC — 

20.00 0.0 

. — 

.60. 00 0.0 


— SOrOO— 0v0™~ 

100.00 0.0 

-l- 26 -r 00 (h-e — 


140.00 0.0 
-IGOrOO—- — O tQ 

120.00 0.0 


-200700 OtO-- 

220.00 0.0 

(hrO — 

■ 2GC. 00 0.0 

:“2S0r00 OtO — 

300.00 0.0 

"~3 2 O';"!) 0 "OtG — 


3JjO. 00 


0.0 


■ 330.00 0.0 

" i^ OO. 00- Ov O --- 


HD-.-0^1O —^-^^-^^-14 

“0. 000 ■“0.3£^i C.OOO 

— “ O-:-0 GO 0 3-S 2- — — -0 V 0 0 G— 

“0.000 “0.389 ■ “0.000 

O-vOO-O -G-.-3-7-G 0 rOOO— — 

-0.000 -0.370 0.000 

— 0-;— 0-Q-G — Q , 3 G 1 — O'vO'OD- 

9.000 “0.348 0.000 

— -0 T 0-60- 0 v33d- 0-.-r0 00 

0.009 “0.309 0. 000, 

-0 - -OOil— -0-.”2-8-2~- 0 -.-0-0 6 

• 0.000 “0.253 “0.000 

“G~r'2'-2-3 6-;-G-9-G 

0.000 “0.194 0.000 

-O-rO-OG “{h-l-GT Gt-OOG 

“0.000 “0.144 0.000 

OtG 00 e-i-l-G-4 OvGOG— — 

0.000 .-0.107 0.000 

OvGOG “-O'rGGG- 0-rGGG — — 

0.000 “0.030 0.000 

0 ; 0 0 0— 0-;-G 7 0 0-.-0 00 
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SATELLITE H1 

RADIUS 

DEPTH 

MAGNETIZATI ON 

INTENSITY 
DECLINATION 
INCLINATION 


= no Km 
= 200 Km 
4 0 K m 


= lO’^ emu 
= 0.0 deg 

= 90.0 deg 


X 



CRATER 


SATELLITE POSITION 


ANOMALY 


X 


y 


AX AY AZ 


0 : 0 0 . O' . 

20.00 0.0 

— 4 0v00 

GO. 00 0.0 

.QQ 

100.00 0.0 

"-1-2OT-O0 -OtO — 

. 140.00 0.0 

-"160.00 -OtO — 

■ 180.00 0.0 

- 200; 00- — OtO — 

220.00 0.0 

—2Wd-r^Q 0 - 0 — 

. 260.00 0.0 

-2SO-rOO“^ OtO-" 

300.00 0.0 
-3 2 0t 00 — — 0-0— 

340.00 0.0 

— 36 0— 0 0 G^H) — 

350.00 0.0 

-400-00 Oi-0- 


0. OOD '0 . 000 0. 7G8 

-0.062 0.000 0.7G7 

“0-1-06 0-.-00 0 e-7 61- 

-0.163 -0.000 0.750 

. -O-rZ-a^ 0-000 0 t-7-3 2~- 

-0.286 0.000 0.702 

-Ot- 5-5-1 —0—06-0 -e-r-G-S-T- 

-0.413 0.000 0.592 

0-t463 0.-000 0.-504- 

-0.491 0.000 0.397 

0 .-4 89 OtOOO — — — Ov-2 8-2- 

-0.457 0.000 0.174 

Ot"^ H) 3 Ot'OO'O 0-j-6“8"7 — 

-0.341 0.000 0.025 

•— 0 -;-2 B-I— O^O 0 0 0-.-0-1-4- 

-0.228 0.000 -0.037 

O-j'1-8 4- O v 0 0 0 • 0 49— 

-0.149 0.000 -0.053 

— O-.-O-O-O — — — * 0- v-6-5-4 — 

-9.098 0.000 -0.052 

Q....Q Q.Q v-04-9- 
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SATELLITE 

RADIUS 

DEPTH 


HT. = ilO Km 

= SOO Km 

= 100 Km' 

MAGNETIZATION 

INTENSITY = 10"^ emu 

DECLINATION = 0.0 deg 

INCLINATION = 0.0 deg 



CRATER 


SATELLITE 

POSITION 


ANOMALY 


K 

X 

AX 

AY 

AZ 

0 » 0 

0.0 • 

- 0.906 

^oTooo 

“ 0.000 

20.00 

0.0 

“ 0.903 

0.000 

“ 0.122 

40.00 

0.0 : 

“ 0.893 

“ 0.000 

“ 0.248 

60 . 00 

0.0 ’ 

“ 0.874 

0 , 000 

“ 0.379 

< 80.00 

0.0 

“ 0.841 

0.000 

“ 0.519 

100.00 

D.O 

“ 0.788 

0.000 

- 0.666 

. 120.00 

0.0 

“ 0. 706 

0.000 

“ 0.818 

140.00 

0.0 

- 0. 586 

0 . 000 

“ 0.965 

160.00 

0.0 

“ 0.422 

0 . 000 

- 1.086 

180.00 

0.0 

- 0,222 

0 . 000 

“ 1.156 

200. 00 

0. 0 

“ 0.012 

0.000 

“ 1.153 

220.00 

0.0 

0.174 

0.000 

“ 1^077 

240.00 

0,0 

0.307 

“ 0. 000 

“ 0.951 

260. 00 

0.0 

0.384 

“ 0.000 

“ 0.807 

280.00 

o . o " 

0.414 

0 . 000 

“ 0. 663 

300 . on 

0.0 

0.415 

“ 0.000 

“ 0.547 

320. 00 

0.0 

0.395 

0.000 

“ 0.446 

340.00 

0.0 

0.368 

“ 0.000 

“ 0.365 

360. 00 

0.0 

0.537 

0.000 

“ 0.300 ■ 

380. 00 

0.0 

< 0.307 

0.000 

“ 0.247 

■ 400 . CO 

0.0 

0.278 

0.000 

“ 0.206 
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SATELLITE HT. = 110 Km 

RADIUS ' = 200Km 

DEPTH = 100 Km 

MAGNETIZATION 

INTENSITY = 10"® emu 

DECLINATION 90.0 deg 

INCLINATION = 0.0 deg 


SATELLITE POSITION 
X y 


0.0 

0.0 

20.00 

0.0 

„ 40 ,. AO 


60.00 

0.0 . 

80.00 


100.^00 

0.0 

120. '00 

0 . 0 

140.00 

0.0 

160.00 

0.0 

180.00 

0.0 

200.00 

0.0 

220.00 

0.0 

240.00 

n.o 

260.00 

0.0 

280.00 

0,0 

300,00 

0.0 

520,00 

0.0 

340.00 

0.0 

360. 00 

0.0 

580.00 

0.0 

400.00 

0.0 


X 



CRATER 

ANOMALY 


AX 

AY 

AZ 


“0.obo„ 

-0.000 

-Q,_aQ.Q..._... 

.-0.906 
-0.905 
~ n n 0 

-0.000 

■ .0.000 
~ fi A A n 


' 0.000 
0_.J}.CA 

-0.896 
- 0 .887. . 

-0.000 
n , nnn 


-0.000 

-Q. OOP 

-0. 873 
~n . 85? 

0.000 
n nnn 


0.000 

^0-,.,i).0.0_._^„ 

-0.823 

O.COO 
- -i3.00Q_ 


0.000 
Q.._noo 

■ -0.733 
-0.071 

0.000 
-0. non 


0.000 

-0.000 

-0. 002 
__-0.3 3 2 

“0.000 
~o - non 


0.000 
o^_Dao 

-0.464 
. -0.407 

0. 000 
fi . non 

i 

-0.000 
0.0 00 

-0.547 

-0.500 

0.000 

“0_^000 


-0.000 
0. non 

. -0.260 
— ril.226 

0.000 
0 . nnn 


0.000 

Q-wQ^OQ 

“0.197 

“Q.175 

0. 000 
n nnn 
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SATELLITE HT. = 110 Km 

RADIUS ■ = 200 Km 

DEPTH = 100 Km 

MAGMETIZATION 

INTENSITY = tO~® emu 
DECLINATION = 0.0 deg 

INCLINATION = 90,0 deg 


X 



CRATER 


SATELLITE 

POSITION 


ANOMALY 


X 

X 

AX 

AY 

AZ 

*- 

0.0 

' ’o.o'"" 

ij . w ij 

- u « 0 0 


20.00 

o.c ■' 

“0.122 

0.000 

1.803 

<j0. 00 ' 

0.0 

-0.24 8 

0 ._.0,0 0,„. „ 


60.00 

0.0 

-0.379 

-0.000 

1.770 

80.00 

0.0 

-0.319 


1....7.2I. 

100.00 

0.0 

-0. G6G 

0.000 

. 1.660 

120.00 

__n.J 

- 0 .^8 1 8. 

0 . ono 

1 .558 

l/}0'. 00 

0.0 

-0.9G5 

0.000 

1.409 

160.00 

0.0 

-1.086 

0,000 

1.206 

IRQ. 00 

b. 0 

-1.156 

0.000 

0.955 

200. 00 

0.0 . 

-1.153 

o.noo 

0.68 3 

220. 00 

o\b 

-1.077 

0.000 

0 . 429 

240.00 

0.0 

-0 J151 

0.000 

0.225 

260.00 

o’. 0 

i -0.807 

0.000 

0.080 

280. 00_ 

n. 0 

-0.668 

0. 000 

-0.012 • 

300 . bo 

0.0 

-0.547 

0.000 

-0 . 066 

5 20 ._o q 

0.0 

-0.446 

0.000 

-0.09,5 

3 40. 0 b 

b. 0 

• -0.365 

0.000 

-0.108 

_ _5 6jl . q q 

0.0 

r0..iqo 

o.noo 

-dJLXL. 

380. bb 

0 . 0 

-0.247 

0.000 

-0.110 

4oq._oq 

q.,q„... 

-q... 2.0,6 

0..QQQ 

_„-0.^LO3. 
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